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We carry out the first complete calculation of kinematic power corrections ~ t/Q 2 and ~ rn 2 /Q 2 
to the helicity amplitudes of deeply-virtual Compton scattering. This result removes an important 
source of uncertainties in the QCD predictions for intermediate momentum transfers Q 2 ~ 1 — 
10 GeV 2 that are accessible in the existing and planned experiments. In particular the finite-t 
corrections are significant and must be taken into account in the data analysis. 
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Deeply Virtual Compton Scattering (DVCS) is the 
simplest process that gives access to generalized par- 
ton distributions (GPDs) and is receiving a lot of atten- 
tion yj, |2| . The existing experimental results come from 
HERMES and Jefferson Lab (Hall A and CLAS) and 
many more measurements are planned after the Jeffer- 
son Lab 12 GeV upgrade and at COMPASS-H at CERN. 
Since the bulk of the existing and expected data is for 
photon virtualities Q 2 < 5 GeV 2 , corrections of the type 
m 2 /Q 2 , t/Q 2 , where m is the target (nucleon) mass and 
t = {p' —p) 2 is the momentum transfer to the target, can 
have significant impact on the data analysis and should 
be taken into account. The finite— t corrections are of 
particular importance if one wants to study the three- 
dimensional picture of the proton in longitudinal and 
transverse plane Q , in which case the ^-dependence has 
to be measured in a sufficiently broad range. 

The necessity of taking into account kinematic power 
corrections to DVCS is widely acknowledged @, I34l3|. 
Early attempts to calculate such corrections by analogy 
to Nachtmann corrections [3] to the structure functions 
in deep-inelastic lepton-nucleon scattering produced re- 
sults that were not gauge invariant and not transla- 
tion invariant with respect to the choice of the posi- 
tions of the electromagnetic currents. The reason is that 
in addition to Nachtmann-type contributions related to 
subtraction of traces in the leading-twist operators one 
must take into account their higher-twist descendants ob- 
tained by adding total derivatives: 0\ ~ d 2 O fJil ... fJin , and 
2 ~ 9' tl O /il ... /lB , where C M1 ... M „ are the usual leading- 
twist operators. The problem arises because matrix el- 
ements of the operator Oi on free quarks vanish [lEj]. 
Thus in order to find its leading- order coefficient function 
in the operator product expansion (OPE) of two electro- 
magnetic currents one is forced to consider either more 
complicated (quark-antiquark-gluon) matrix elements, or 
stay with the quark- antiquark ones but go over to the 
next-to-leading order in a s . In both cases the real diffi- 
culty is not the calculation of the relevant Feynman dia- 
grams, but the necessity to separate the contribution of 
interest from the "genuine" quark-gluon twist-four oper- 



ators. 

This problem was solved in Refs. 0, [TtJ using confor- 
mal symmetry which implies that coefficient functions of 
"kinematic" and "genuine" twist-four operators are mu- 
tually orthogonal with a proper weight function [Lsl. [l9| . 
Using this approach we have calculated in Ref. [19J the 
finite-i and target-mass corrections to DVCS for the 
study case of a scalar target. We verified gauge- and 
translation-invariance and, most importantly, found that 
the structure of kinematic corrections proves to be consis- 
tent with collinear factorization. In this letter we present 
our final results for the helicity amplitudes of DVCS to 
the 1/Q 2 accuracy for the physically interesting case of 
the spin-1/2 (nucleon) target. This result removes one 
important source of uncertainties in the QCD predictions 
for intermediate photon virtualities that are accessible in 
the existing and planned experiments. 

The DVCS amplitude j*(q) + N(p) -> -y(q') + N(p') is 
defined by the matrix element of the time-ordered prod- 
uct of two electromagnetic currents, sandwiched between 
the nucleon states 

ijd^xj d 4 ye-^ + ^^(p'|T{^ m (x)jr(y)}|P> = 

= {2TT) 4 S( P + q-p'-q')A ftu (q,q',p). (1) 

Introducing the photon polarization vectors (e ± '°g) = 0, 
{e q') = one can write A^ in terms of helicity ampli- 
tudes 



A — F + F ~A ++ 

+ 4e+A°- 



£fj, A 



e°eZA 0+ 



-e , e, 



lAf. (2) 



The last term ~ is of no interest as it does not con- 
tribute to any observable. 

The helicity-conserving amplitudes are the leading 
ones in the scaling limit, A ±±: ~ O(Q ), and the hclicity- 
fiip amplitudes are power-suppressed: A a± ~ 0(Q^ 1 ), 
A ±T ~ (D(Q~ 2 ). Thus in order to calculate physical 
observables to the 1/Q 2 accuracy one has to take into 
account 1/Q 2 corrections to A ++ and A~~, whereas for 



2 



the helicity-flip amplitudes the leading power accuracy is 
sufficient. 

The definition of helicity amplitudes depends on a ref- 
erence frame. We use the photon momenta, q and q', to 
define a longitudinal plane spanned by the two light-like 
vectors 

n = q', h= -q+(l-T)q', (3) 

where t = t/(Q 2 + t), Q 2 = —q 2 . For this choice the 
momentum transfer to the target A = p' — p = q — q', 
t = A is purely longitudinal and the target (proton) 
momenta have a nonzero transverse component 

\Px\ 2 = -m 2 - „ O(Q ) , (4) 



GPDs Hi, Ei, Hi, Ei. For the GPD definitions we follow 
Ref. [J. 

The calculation is similar to the case of the scalar tar- 
get considered in Ref. [l9[ so that in this letter we only 
present the final expressions. Note that the electromag- 
netic gauge invariance is guaranteed to twist-four accu- 
racy already on the operator level and is embedded in 
the definition of helicity amplitudes. The translation in- 
variance (independence on the shift of the positions of the 
electromagnetic currents in Eq. ([1]): x — > x + S, y — > y+S) 
is nontrivial and provides a strong check of the calcula- 
tion, see Ref. [l9j]. The results can conveniently be writ- 
ten in terms of the vector and axial- vector bispinors 

V** = u(p')Yu{p) , a? = u{p')Yl5u{p) . (7) 



where P — (p + p')/2 and the skewedness parameter is 
defined as £ = -(A • q')/(2(P ■ q')) . 

The condition |PjJ 2 > translates to the lower bound 

\t\>\t min \ = 4m 2 e/(i-e),ci 0. 

We choose the polarization vectors as follows [19| 



We define v ± — (u-e d 



(o- e ±), P± = (P-£ d 



Although P+ = P^ = \P±\/y/2 we prefer this notation 
to keep trace of the polarization vectors. We also use a 
shorthand notation 



(q» -%q 2 /(qq')) /v 7 ^? 2 - 



(5) 



X*(x,£,t) =H«(x,£,t)+E 9 (x,t;,t) 



(8) 



and rewrite the helicity-conserving amplitudes in terms 
of the vector- and axial- vector invariant functions as 



where P± = g^P", P± = ej^P" and 

atu =9»u - {%<( v + q'^Iu)/ (qq') + q'^ql q 2 l (qq') 2 , 

4u =(^q a q' P /(qq') ■ (6) 

Each helicity amplitude involves the sum over quark 
flavors, A = Yl e qAq, where e q is the quark electromag- 
netic charge, and is written in terms of the leading- twist 



\{A + q + +A- q . 
1 -{A+ + -A--) = 



(vP) 
2m 2 
(aA) 
Am 2 



M 

(qq') 
M 
(qq') 



v 2 J 



(9) 
(10) 



The following expressions for V^, A^ present our main 
result: 



V? = 




t 
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VI = 
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A? = 
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w 




t 1 - 




A| = 





E" ® Cn" 



Q 2 



E q ® Cf 



a (t + 2 \p ± n%)edtE« ® c 2 - + ^feo^x* ® c 2 - 



Q 2 



4 

Q~ 2 



(11a) 
(lib) 



2 ft 
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A ? = 1 - ^)iE q ® c+ + —ze« ® c+ - — - + 2|p ± i 2 r^k 2 ^e-B 9 ® c. 



£H q ®C+ + —£H q ®C+ 
Us 



y 2 



8m z 
£P«® C+ 



£ 2 <%P 9 ® C 2 + , (11c) 
(lid) 



r 



In addition, for the helicity-flip amplitudes we obtain 



Q 2 



± [af - 4Pf ® C? 



Pi 



: («P) 



e 2 <% p« ® cr 



2m z 



(12) 



and 



8P? 



vf - 2P 



x Q 2 



e 2 a 5 ra 5 x«® [zcr 



± (a 

+ p 



2Pf 



±(^) ,3 



edAed^H q ®c+ 



(13) 
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In all cases d$ = d/d^ and the notation F®C stands for 
the convolution of a GPD F with a coefficient function 
C: 

F®C = J dxF(x,£,t)C(x,£) . 

The coefficient functions C^(x,£) are given by the fol- 
lowing expressions: 



1 



1 



± 



■In 



le £ — x 
2£ 



Li 2 



2i 



iej ± (x O —x) , 
-Lia(l) 



±(a:o-a;)^ + -Cf(a;,0- 



(14) 



Note that Cq have simple poles at i = ±£ whereas 
C^ 2 have a milder (logarithmic) singularity at the same 
points. This ensures that the kinematic power correc- 
tions are factorizable, at least to the leading order in a s . 
In the DVCS kinematics 

{vq') ~ (aq') = 0(Q 2 ) , (vP) ~ (oA) = O(Q ) (15) 



so that the helicity-conserving amplitudes (|lla|) - (|lld[) 
include leading contributions O(l/Q ) and the correc- 
tions C(l/(5 2 ), whereas all terms in Eqs. (JT2J) and (JTSJ) 
are of the order 0(1/ Q) and 0(1/Q 2 ), respectively, as 
expected. 

Our results for the DVCS helicity amplitudes have a 
similar structure to those in Ref. [19( for the scalar tar- 
get [2(| ■ The main difference is the appearance of a large 
target mass correction to the contribution of the GPD 
E q (E q ) that involves X q (H q ), cf. the last term in the 
first line of Eq. (TTTal) (Eq. (fTTc|> ). 

It has become customary to parametrize the DVCS 
amplitude @ in terms of the so-called Compton form 
factors (CFFs) % and £ @ 



2m 



In our notation 

w = Vi-v 2 /e, 



(ua^q' A v u) £ 



£ = -Vi 



(16) 
(17) 



where \ 

A detailed study of the numerical impact of the kine- 
matic corrections on different DVCS observables goes be- 
yond the tasks of this letter. For orientation, we have 
calculated the corrections to the imaginary parts of the 
CFFs H and £ which involve the GPDs in the DGLAP 
region only To this end we use the GPD model of 
Refs. ElLlil: 



E q 



(x, & t) = dp da 6(x-j3- £a)f q (p, a, t) , 
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FIG. 1: The coefficients (£,t) (upper panel) and e^(£,t) 
(lower panel) for different values of t (in GeV 2 ). 



where the double distribution f q is written as 



f q (l3,a,t) = h(/3,a)\l3\- a ' t q((3) 



Here q(/3) is the MRST2002 NNLO valence u- and d- 
quark distribution [23[ ) and the profile function h is given 
by the following expression: 



« 2 )/(l-|/3|) 3 



where k u ~ 1.7 and fc^ ~ —2.0 are the anomalous 

magnetic moments, r\ u ~ 1.7 and r/^ ~ 0.57 |22l |. and 

A, 

ratios 



q — Jdf3(l — (3) nq q(f3) . We consider the following 



ImJ r -ImJ rLO t r _ m 2 



= 7^ c f(^) + 7^ c ™(^) 



ImJ" LO 
where T = {%,£}, 

imH LO =tt £ e 2 [ffn-e, e, t) - # 9 (£, £, *)] 



(18) 



(19) 



and similar for Im£ L °. The coefficients c[(£,t) depend 
on t because of the non-factorizable ^-dependence of the 
GPDs through the Regge trajectory. 

In Fig. Q] we show c^(£,i) and c^(£, t) as a function 
of £ for several ^-values: -t = 0,0.5, 1.0, 1.5 GeV 2 . The 



4 



-0.7 




'0.0 0.1 0.2 0.3 0.4 0.5 




0.0 0.1 0.2 0.3 0.4 0.5 

S 

FIG. 2: The coefficients cf(£,t) (upper panel) and c^(£, t) 
(lower panel) for different values of t (in GeV 2 ). 



same is shown in Fig.[5]for cf (£, i) and c~j(£, t). One sees 
that the corrections to £ are in general larger than for 
the % form factor; in particular £ receives a relatively 
large proton mass correction. 

Finally, note that the finite-i correction depends on 
the definition of the skewedness parameter £, which is 
not unique. If one defines £ through the Bjorken xb 
parameter, £g = £s/(2 — xb), which seems to be natural 
from the experimental point of view (the relation of "our" 
£ to £ B is given in Eq. (125) in Ref. [II]), c^' f (£, t) change 
accordingly, but in general do not become smaller. 

To summarize, in this work we have calculated, for 
the first time, the kinematic power corrections ~ t/Q 2 
and ~ m 2 /Q 2 to the helicity amplitudes of deeply virtual 
Compton scattering. These corrections are important for 
intermediate momentum transfers Q 2 ~ 1 — 10 GeV 2 that 
are accessible in the existing and planned experiments, 
and have to be taken into account in the data analysis. 
In particular the finite— t corrections are indispensable if 
one aims to study "holographic" images of the proton 
in the transverse plane, in which case the i-dependence 
must be measured in a broad range. 
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